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Abstract 

Some aspects of the theory of cosmological perturbations from cosmic strings 
and other topological defects are outlined, with particular reference to a sim- 
ple example: a spatially flat CDM-dominated universe. The conserved energy- 
momentum pseudo-tensor is introduced, and the equation for the density pertur- 
bation derived from it. It is shown how the scaling hypothesis for defect evolution 
results in a Harrison-Zel'dovich spectrum for wavelengths well inside the horizon. 



*From: Proceedings of 'Trends in Astroparticle Physics', Stockholm, Sweden 22-25 September 1994, 
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One of the most exciting problems in cosmology today is to explain the origin of the 
primordial fluctuations, whose imprint we have seen in the temperature of the cosmic 
microwave background 0. It is generally thought that these fluctuations owe their 
origin to processes happening very early on in the history of the Universe, perhaps at 
about 10~ 36 seconds after the Big Bang. In the last 15 years or so, two classes of 
theory have emerged to explain how this came about: those based on inflation and those 
on ordering dynamics after a cosmological phase transition (sometimes loosely called 
topological defects). 

Inflationary fluctuations (see e.g. |2|) arise from the amplification of quantum fluctua- 
tions in a weakly-coupled scalar field by the strong gravitational field of an early phase of 
accelerated expansion. They have the statistical properties of a Gaussian random field, 
for which the calculations of correlation functions (such as the power spectrum of density 
fluctuations) are relatively easy. 

Cosmic strings || and other defects [|J, on the other hand, are certainly not Gaussian. 
They are formed in an early phase transition, in an initial distribution more akin to a 
Poisson process, although one of the major problems of the string scenario is to under- 
stand the statistical properties of the evolving network. In this talk I shall show how 
an understanding of string statistics can be converted into knowledge of the statistics 
of the primordial fluctuations, using the general relativistic theory of cosmological per- 
turbations from "stiff" sources developed by Veeraraghavan and Stebbins ||. To bring 
out the essentials I shall consider a very simple case: an Q = 1 universe dominated by 
cold dark matter (CDM), and study the density perturbation only. Much of what I say 
applies to textures || as well. 

Firstly, I outline notation and conventions. For the space-time metric I write 

g»u(v, x ) = a2 (v)(v^ + h^uiv, aO), (1) 

so that r) is conformal time, x are the comoving coordinates, r)^ = diag(— 1, 1, 1, 1, ) 
is the Minkowksi metric, and is the (scaled) metric perturbation. I shall use the 
synchronous gauge ho^ = 0. Greek indices from near the middle of the alphabet run from 
to 3, while Roman ones are purely spatial indices. Natural units are used throughout. 

Strings are stiff sources: that is, their energy-momentum tensor is separately con- 
served, so that 

e 00 + -8+ + diQoi = 0, (2) 
a 

Qoi + -©oi + djBij = 0. (3) 

CL 

These equations follow from the equations of motion of the sources, which are unaf- 
fected by small variations in the background gravitational field. For strings, can 
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be expressed in terms of the linear mass density /x and the coordinates of the string 
world-sheet, X^(a, rf): 

9^(ry, x)=njdo {X^X V - X"X v )5®{x - X(a, r,)), (4) 

where X^ = dX^/da and = dX^/drj. 

The evolution of the perturbations to the metric and to the energy-momentum tensor 
of the matter content of the Universe are then described by the linearised Einstein equa- 
tions |7| , together with the covariant energy- momentum conservation of the matter. This 
evolution can be neatly expressed as the ordinary conservation of an energy-momentum 
pseudo-tensor t^ v [[5]. For CDM plus strings (or other stiff sources) it takes the form 

r 00 = aV-^A + Qoo 
07rG a 



r° 4 = -0 



o? 

T " = -dcsfa - !**«) + e - (5) 

where p is the background CDM density, 6 is the fractional perturbation to it, and 
is the traceless part of hij. The conservation equation d^r^ = reduces in this case, 
after the use of the equations of motion of the defects (S) and the constraint equation 
h = -26, to 

+ - * - 5 (-V * = 47TG9+, (6) 

where 6 + = O o + 6n- This is a second order inhomogeneous linear ODE, and is easily 
solved in terms of two Green's functions, 

GM) = i(3^ + 2^), (7) 
5 \ rf rf ) 

given a source function Q + (rj,x) and initial conditions 6i and 6i at time rji. The initial 
conditions are not, arbitrary, for they must have arisen from a causal evolution, satisfying 
pseudo-energy-momentum conservation. On super-horizon scales at rji, the total pseudo- 
energy vanishes, or 

a 2 lPl 6i{k) + -r^-K + 6 00 = 0. (9) 
47rG a 

This equation expresses the fact that energy is locally conserved when making strings (or 
other defects). The initial energy in the strings is compensated by that in the matter. 
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The general solution of equation (|6]) can be split into two parts, as is usual with 
inhomogeneous differential equations, 



^{n) = Gi(r],rii)6i + Q2(r] 1 rii)6i 1 (10) 
5 s (rj) = fdrfg^fjDQ^rf). (11) 



These were termed the initial and subsequent compensation by Veeraraghavan and Steb- 
bins 0. The initial compensation soon settles into a purely growing mode, with 
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^oovw (12) 

where I have used @ in the second equality. We see that the initial compensation is 
entirely determined by the initial defect density. 

This raises an important question: how does the Harrison-Zel'dovich spectrum (|5| 2 ) ~ 
rfk arise? This spectrum is special because it is scale invariant: the r.m.s. mass fluc- 
tuations in a horizon-sized volume are constant. To demonstrate how a scale-invariant 
spectrum emerges from the assumption of scaling in the network, we must first show that 
the initial compensation is cancelled by a piece of the subsequent compensation. After 
all, we should not expect the power spectrum to depend on the initial conditions, for 
they may well produce white noise (k°), having been generated by a causal process inside 
the horizon at rji. Using the equation of motion of the defect network (|3|), we find 




2 J m r/ 3 

+ ToJj^^-lf) e M- (13) 

The first term on the first line cancels the initial compensation, the second is a transient, 
and we are left with the second line as the true source of the density perturbations in 
the dark matter component. Thus it is extremely important to get the initial conditions 
right in numerical simulations, otherwise one may well find spurious white noise in the 
final power spectrum, which will dominate any Harrison-Zel'dovich component at late 
times |J. 

To calculate the power spectrum, we need the density perturbations well inside the 
horizon, for which kr] ^> 1. All that survives (provided Qq^t], k) drops off fast enough) 
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is the growing mode, which is the piece proportional to rf. Thus 

(m\ 2 ) - (^) 2 v%kj £dvidv2 (e w (i fc) fc)es i (i b) fc)>. (u) 

The angle brackets denote an average over an ensemble of realisations of the defects. We 
can deduce the form of the right hand side, using the scaling hypothesis for defects. 

Firstly, we do not expect the momenta Qoi(rj,k) to be correlated at two different 
times T]i, T] 2 if 1 772 — Vi\ ^ k~ x . Let us denote the correlation time by Tk^ 1 . Secondly, 
statistical isotropy in the defect network means that (©oifa^cyfa)) = ^%(l©oi(^)| 2 )- 
Thus, 

1 / AttO\ 2 rv 

mv)\ 2 ) - 3 (^-) v 4 rk jf d % (leoifa.fe)! 2 ). (is) 

Thirdly, the scale invariance tells us that the r.m.s. momentum density fluctuations are 
constant on the horizon scale: 

^^(l©o,(^,fc)| 2 )oc^ 4 . (16) 

Therefore we may write 

(\Q 0l (r ] ,k)\ 2 ) = V- 1 r ] - 1 ^S^\kr ] ), (17) 

where S^ 0l \z) is a dimensionless scaling function associated with the momentum density, 
and V is the normalisation volume. Finally, we have 

<l%)| 2 > - V-> 1 - (±^) Vt/c/J^S^). (18) 

Thus, provided the momentum scaling function is well-behaved, and the integral con- 
verges, we recover the Harrison-Zel'dovich spectrum. There is perhaps a logarithmic 
correction, if the high frequency limit of the source power spectrum is white noise, as 
would be supplied by a population of small string loops. Given the expression for Q 0i , 
easily calculable from (|j), we should expect to be proportional to the mean square 
string velocity v 2 and the average number of Hubble lengths of string per Hubble volume 

7 2 - 

A similar procedure was adopted by Albrecht and Stebbins to derive an analytic 
approximation to the power spectrum for strings. However, they did not explicitly remove 
the initial compensation, and simulated it by a window function. One then finds that 
the power spectrum depends on (|6 + | 2 ), which goes as v A . Nevertheless, the functional 
form of their approximation that they give is reasonable, even if the details are not quite 
right. 
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In general terms, we learn from the above that the appropriate objects for the nu- 
merical study of the density perturbation power spectrum are the two-time correlation 
functions (©oo^i)©^^)) and (0oj(^i)@oi( ? ?2))- Getting good statistics is rather hard 



with fully realistic string simulations ||, but simulations in flat space |10[ run very 



much faster, and are thus ideal for numerical experiments on the form of the correlation 



functions, and the calculations of the observable quantities |l T|. It is to be hoped that 
they will also aid us in understanding the coarse-grained features of string evolution and 
in constructing analytic models p^| , which attempt to derive dynamical equations for 
quantities such as 7 and v. 
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